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We report on a search for pair production of first-generation scalar leptoquarks (LQ) in pp colli- 
sions at y / s=1.96 TeV using an integrated luminosity of 203 pb~ x collected at the Fermilab Tevatron 
collider by the CDF experiment. We observe no evidence for LQ production in the topologies arising 
from LQLQ — * eqeq and LQLQ — > eqvq, and derive 95% C.L. upper limits on the LQ production 
cross section. The results are combined with those obtained from a separately reported CDF search 
in the topology arising from LQLQ — > vqvq and 95% C.L. lower limits on the LQ mass as a function 
of (3 = BR(LQ -> eq) are derived. The limits are 236, 205 and 145 GeV/c 2 for /3 = 1, /3 = 0.5 and 
— 0.1, respectively. 

PACS numbers: 14.80.-j, 13.85.Rm 



The remarkable symmetry between quarks and leptons 
in the Standard Model (SM) suggests that some more 
fundamental theory may exist, which allows for new in- 
teractions between them. Such interactions are mediated 
by a new type of particle, the leptoquark (LQ)Q] and are 
predicted in many extensions of the SM, e.g. grand uni- 
fication, technicolor, and supersymmetry with R-parity 
violation^. A LQ carries both lepton and baryon num- 
ber, is a color triplet boson with spin or 1, and has 
fractional charge. Usually it is assumed that LQs couple 
to fermions of the same generation to accomodate exper- 
imental constraints on flavor changing neutral currents 
and hclicity suppressed decays. 

Previous experimental limits on LQ production are 
summarized in 3]. The HI and ZEUS experiments at the 
e ± p collider HERA published 4lj lower limits on the mass 
of a first generation LQ that depend on the unknown 
LQ I — q Yukawa coupling A. At the LEP collider, pair 
production of LQs can occur in e + e~ collisions via a vir- 
tual 7 or Z boson in the s— channel and lower limits have 
been presented in [jj. At the Fermilab Tevatron Q 
LQ would be predominantly pair produced through qq 
annihilation and gg fusion. Since the production is me- 
diated via the strong interaction it is independent of A, 
in contrast to the searches at e-p machines. The coupling 
strength to gluons is determined by color charges of the 
particles, and is model-independent in the case of scalar 
LQs. The production of vector LQ pairs depends on ad- 
ditional assumptions on LQ coupling to gluons and its 
cross section is typically larger than the cross section for 



scalar LQs production. Since the acceptance for vector 
and scalar LQ detection is similar, limits on the vector 
LQ mass will be more stringent. 

In this Letter, we focus on a search for first-generation 
scalar LQ pairs produced in pp collisions at ^=1.96 
TeV. A search for scalar LQ pairs decaying into vvqq, 
resulting in jets and missing transverse energy topology 
has been presented in ||. Here we study alternative fi- 
nal state signatures, with LQs decaying into eejj and a 
final state consisting of two electrons and two jets and 
LQs decaying in evjj and the final state consisting of an 
electron, two jets, and missing transverse energy . The 
results for three channels are combined and presented as 
a function of /?, the LQ branching fraction into an elec- 
tron and a quark. 

CDF is a general-purpose detector built to study the 
physics of pp collisions at the Tevatron accelerator at Fer- 
milab and it is described in detail in The data used in 
the analysis were collected during the 2002-2003 Tevatron 
Run II. The integrated luminosity for this data sample 
is 203 ± 12.2 pb _1 . Events are selected if they pass the 
high Et electron trigger, requiring one electromagnetic 
trigger tower to be above threshold and a set of iden- 
tification cuts on the electromagnetic cluster, track and 
shower profile. The efficiency of the trigger combinations 
used in the eejj and evjj analyses has been measured 
using Z -> ee data[H E3 and it is - 100%. Elec- 
trons are reconstructed offline as calorimeter electromag- 
netic clusters matching a track in the central-tracking 
system (central electrons, \rj\ < l.O^H) or as calorime- 
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ter electromagnetic clusters only in the forward region 
(1 < \v\ < 3). Electromagnetic clusters are identified 
by the characteristics of their energy deposition in the 
calorimeter: cuts are applied on the fraction of the en- 
ergy in the electromagnetic calorimeter and the isolation 
of the cluster. The identification efficiency for a pair of 
central electrons is 92.4% ± 0.4 and for a pair of central- 
forward electrons is 80% ± 0.4. The coordinate of the 
lepton (also assumed to be the event coordinate ) along 
the beamline must fall within 60 cm of the center of the 
detector ( z vertex cut) to ensure a good energy measure- 
ment in the calorimeter. This cut has an efficiency of 
95% ± O.l(stat) ± Q.5(sys), and it was determined from 
studies with minimum bias events. The efficiencies of 
the identification cuts, the trigger selection and the ver- 
tex cut, measured using Z — ► ee data were taken into 
account when evaluating the signal acceptance and back- 
ground estimate. Jets are reconstructed using a cone of 
fixed radius R = y (Arj 2 + Ac/) 2 ) = 0.7 and required to 
have 1 77 1 < 2.0. Jets have been calibrated as a function 
of r\ and Et and their energy is corrected to the parton 
level|l2|. Neutrinos produce missing transverse energy, 
I£t, which is measured by balancing the calorimeter en- 
ergy in the transverse plane. 

In the analyses we describe here, the signal selection 
criteria are set according to the kinematic distribution 
(e.g. Et of the electrons and Et of the jets) of de- 
cay products determined from Monte Carlo (MC) stud- 
ies, optimized to eliminate background with a minimal 
loss of signal events [1^ Hif . In the dielectron and jets 
topology, we select events with two reconstructed iso- 
lated electrons with Et > 25 GeV. At least one electron 
is required to be central, while the other can be central 
or forward. Events are further selected if there are at 
least two jets with Et > 30 and 15 GeV. The dataset 
selected above is dominated by QCD production of Z 
bosons in association with jets and tt production where 
both the W's from top decay into an electron and neu- 
trino. To reduce these backgrounds the following cuts 
are applied: i) veto of events whose reconstructed dilep- 
ton mass falls in the window 76 < m ee < 110 GeV/c 2 
to remove the most of the Z + jets contribution, ii) 
E T (ji) + E rih) > 85 GeV and E T (e x ) + E T (e 2 ) > 85 
GeV, iii) y/(Er{ji) + ErWfi 2 + (E T (ei) + E T (e 2 )) 2 > 
200 GeV to remove the remaining Z + jets and top 
contributions. We studied the properties of the physics 
backgrounds by generating the process Z + 2 jets with 
ALPGEN0 + HERWIG (to perform parton show- 
ering) and tt with PYTHIA then passing them 
through a complete simulation of the CDF II detec- 
tor based on GEANT0] and full event reconstruction. 
Other backgrounds from bb, Z — > tt, WW are negligible 
due to the electron isolation and large electron and jet 
transverse energy requirements. To normalize the num- 
ber of simulated events to data we used the theoreti- 
cal cross sections for tt from [2(| and for j/Z — > ee + 



TABLE I: Efficiencies after all cuts with total error (statis- 
tical and systematic) and 95% C.L. upper limits on the pro- 
duction cross section x branching fraction Br, as a function 
of Mlq, for the two channels. 

Mlq (GeV/c 2 ) eejj evjj 



e (%) axBr(pb) e (%) axBr(pb) 



100 


7 ± 0.5 


1.11 


2 ± 0.26 


5.71 


140 


12 ± 0.5 


0.25 


8 ± 0.7 


0.69 


160 


21 ± 0.8 


0.14 


8 ± 0.7 


0.65 


200 


32 ± 1.2 


0.09 


16 ± 1.3 


0.37 


220 


35 ± 2.0 


0.08 


19 ± 1.5 


0.24 


240 


38 ± 2.0 


0.07 


20 ± 1.6 


0.23 


260 


40 ± 2.0 


0.07 


22 ± 1.7 


0.22 



2 jets from |2lJ. The expected number of Z + 2 jets 
events is 1.9 ± 0.4. The expected number of ti events is 
0.35 ±0.06 events. The background arising from multijet 
events where a jet is mismeasured as an electron (fake) is 
calculated using data, for both this analysis and the one 
that follows. The method used relies on the assumption 
that the fake electron produced by a jet will be accom- 
panied by other particles produced by the fragmentation 
of the jet; thus the isolation fraction of the fake electron 
will generally be larger than the one corresponding to a 
real electron. The isolation fraction is defined here as: 
(E T one - E T luster ) /E T luster where E T one is the sum of the 
electromagnetic and hadronic transverse energies mea- 
sured in all towers in a radius R — v/ (A<fi 2 + Arj 2 ) = 0.4 
around the electron and E T luster is the transverse electro- 
magnetic energy of the electron. The phase space corre- 
sponding to the two electron isolation fractions {eejj) or 
to one electron isolation fraction and the Tf*T (evjj) is di- 
vided in different regions. We assume that there is no cor- 
relation between the isolation of the two electrons (eejj) 
and the isolation of the electron and TpT (evjj). In the 
region where both electrons have large isolation fraction 
(eejj), or where the J^t is small and the isolation frac- 
tion of the electron is large (evjj) the LQ contribution 
is expected to be negligible. We call these background- 
dominated regions. With these assumptions from the ra- 
tio of the number of events in the background-dominated 
regions we can extrapolate the contribution in the signal 
region. We estimate O+q' 7 fake events in the central- 
central category and 4.0 ± 2.0 in the central-forward cat- 
egory. The final background estimate is 6.2 ± 2.2 events. 
We checked the prediction of our background sources 
with data in a control region defined by requiring two 
electrons with E T > 25 GeV, 2 jets with E T > 30 GeV 
and 66 < m ee < 110 GeV/c 2 . We observe 107 events in 
agreement with 113 ± 15 predicted from SM processes. 

The efficiency to detect our signal was obtained from 
MC simulated LQ (PYTHIA) to account for kinematical 
and geometrical acceptances and it is reported in Table 
I for various LQ mass values. The following systematic 
uncertainties are considered when calculating signal ac- 
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ceptance and background predictions: luminosity (6%), 
choice of parton distribution functions (2.1%), statisti- 
cal uncertainty of MC (< 1%), jet energy scale (< 1%), 
statistics of Z — > e + e _ sample (0.8%) and z vertex cut 
(0.5%). After all selection cuts, 4 events are left in the 
eejj channel data. 

In the search in the electron and neutrino plus two 
jets topology, we select events with one reconstructed 
isolated electron with Et > 25 GeV. The electron is re- 
quired to be central (\r)\ < 1.0). We veto events with 
a second central or forward electron to be orthogonal 
to the previous analysis. We then select events where 
there is a large missing transverse energy, I£t > 60 GeV 
and at least two jets with Ey > 30 GeV in the range 
1 77 1 < 2. This time the selected dataset is dominated by 
QCD production of W bosons in association with jets 
and top quark pairs, where either both the W's from the 
top pair decay into Iv and one lepton is not identified, or 
one of the Ws decays leptonically and the other hadron- 
ically. A small source of background is represented by 
Z + 2 jets, where one of the electrons is not identified. 
To reduce these backgrounds the following cuts are ap- 
plied: i)A(j)(j^T — jet) > 10° to veto events where the 
transverse missing energy is mis-measured due to a jet- 
pointing to a non instrumented region of the calorimeter 
, ii)£rOi) + Et{]2) > 80 GeV, hi) transverse mass of 
electron- neutrino system, Mt{zv) > 120 GeV/c 2 to re- 
duce the W + 2 jets contribution. We studied the prop- 
erties of the W + jets, ti and Z + 2 jets backgrounds 
using MC simulated events (ALPGEN + HERWIG and 
PYTHIA). The background from W -> v t t + 2 jets 
(ALPGEN+HERWIG) is negligible after the final mass 
window cut (see below), as as is the background com- 
ing from misidentified leptons and false I^t- Our final 
cut consists in selecting events falling in a mass window 
defined around the LQ mass in the following way. We cal- 
culate the invariant mass of the electron-jet system and 
the transverse mass of the neutrino-jet system. Given 
the decay of the two LQs, there are two possible mass 
combinations for the electron and the neutrino with the 
two leading jets. We choose the combinations that mini- 
mize the difference between the electron-jet mass and the 
neutrino-jet transverse mass. We fit the peak of the e- 
jet distribution with a Gaussian, to obtain an estimate of 
the spread of the distribution in the signal region (er e ), as 
well as the v — jet transverse mass distribution, to obtain 
(T, y . In the kinematic plane of m(e — jet) vs niTiv — jet) 
we define the sides of rectangular boxes centered around 
various nominal LQ mass as 3 x <7 e .„. For each LQ mass, 
events are accepted if they fall inside the rectangular box. 
The overall selection efficiency for various LQ masses is 
given in Table I. We checked the simulation prediction of 
our background sources with data in a control region de- 
fined by requiring one electron with Et > 25 GeV, Ifir > 
35 GeV and 2 jets with E T > 30 GeV. We observe 536 
events in agreement with 503 ± 22 predicted from SM 



TABLE II: Final number of events surviving all cuts in the 
electron, missing energy and jets topology, compared with 
background expectations, as a function of the LQ mass (in 
GeV/c 2 ). Errors include statistical and systematic uncer- 
tainty. 

Mass W + 2 jets top Z + 2 jets Lotal Data 
"^120 1.5 ± 0.9 3.3 ± 0.5 0.06 ± 0.01 4.9 ± 1.0 6~ 
140 1.5 ± 0.9 3.1 ± 0.6 0.08 ± 0.02 4.7 ± 1.0 4 
160 2.5 ± 1.1 2.8 ± 0.6 0.08 ± 0.02 5.4 ± 1.2 4 
180 2.5 ± 1.1 2.4 ± 0.5 0.08 ± 0.02 5.0 ± 1.2 4 
200 2.5 ± 1.1 2.0 ± 0.5 0.07 ± 0.02 4.6 ± 1.2 4 
220 2.0 ± 1.0 1.6 ± 0.3 0.06 ± 0.02 3.7 ± 1.1 2 
240 2.0 ± 1.0 1.1 ± 0.3 0.06 ± 0.02 3.1 ± 1.0 2 
260 1.5 ± 1.0 0.8 ± 0.3 0.04 ± 0.02 2.4 ± 0.9 2 



processes. 

The efficiency to detect our signal was obtained from 
MC simulated LQ data (PYTHIA). The following sys- 
tematic uncertainties are considered when calculating 
signal acceptance and background predictions: luminos- 
ity (6%), choice of the parton distribution functions 
(2.1%), statistics of MC (< 1.0%), jet energy scale (< 
1%), electron identification (0.6%), z verte x cut (0.5%), 
initial and final state radiation (1.7%). The number of 
events in each mass region, after all selection cuts, com- 
pared with the background expectations is reported in 
Table II. 

In the analyses described above the number of events 
passing the selection cuts is consistent with the expected 
number of background events. The conclusion of the two 
searches is that there is no LQ signal: hence we derive an 
upper limit on the LQ production cross section at 95% 
confidence level. We use a Bayesian approach [22^ with a 
flat prior for the signal cross section and Gaussian priors 
for acceptance and background uncertainties. The cross 
section limits are tabulated in Table I and the mass limits 
are tabulated in Table III. To compare our experimental 
results with the theoretical expectation, we use the next- 
to-leading order (NLO) cross-section for scalar LQ pah- 
production from [13 with CTEQ4 PDF0- The the- 
oretical uncertainties correspond to the variations from 
Mlq/2 to 2M^q of the renormalization scale fi used in 
the calculation. To set a limit on the LQ mass we com- 
pare our 95% CL upper experimental limit to the theoret- 
ical cross section for /1 = 2Mlq, which is conservative as 
it corresponds to the lower value of the theoretical cross 
section. We find lower limits on M(LQ) at 235 GeV/c 2 
(/3 = 1) and 176 GeV/c 2 (J3 = 0.5). To obtain the best 
limit, we have combined the results from the two decay 
channels just described with the result of a search for 
LQ in the case where the LQ pair decays to neutrino 
and quark with branching ratio Br(LQ — ► vq) = 1.00. 
The individual channel analyses are in fact optimized for 
fixed values of (3 (1,0.5,0) while in the combined analysis, 
due to the contributions of the different decay channels, 
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TABLE III: 95% C.L. lower limits on the first generation 
scalar LQ mass (in GeV/c 2 ), as a function of p. The limit 
from CDF0 (eejj) Run I (~ 120pfe _1 ) is also given. 
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FIG. 1: LQ mass exclusion regions at 95% C.L. as a function 
of Br(LQ^ eq). 

the signal acceptance can be naturally expressed as a 
function of (3. As for the treatment of uncertainties, the 
searches in the eejj and evjj channel use common crite- 
ria and sometimes apply the same kind of requirements 
so the uncertainties in the acceptances have been consid- 
ered correlated. When calculating the limit combination 
including the vvjj channel the uncertainties have been 
considered uncorrelated. For each (3 value a 95% C.L. 
upper limit on the expected number of events is returned 
for each mass, and by comparing this to the theoretical 
expectation, lower limits on the LQ mass are set. The 
combined limit as a function of (3 is shown in Figure 1, to- 
gether with the individual channel limits. The combined 
mass limits are also tabulated in Table III. 

In conclusion, we have performed a search for pair pro- 
duction of first generation scalar LQs using 203 pb _1 
of proton-antiproton collision data recorded by the CDF 
experiment during Run II of the Tevatron. The results 
from all the final state signatures are combined and no 
evidence of LQs production is observed. Assuming that 



a scalar LQ decays to electron and quark with variable 
branching ratio f3 we exclude LQs with masses below 236 
GeV/c 2 for j3 = 1, 205 GeV/c 2 for (i = 0.5 and 145 
GeV/c 2 for (3 = 0.1. 
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